System-size dependence of A and H production in nucleus— nucleus collisions at 40A 
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Results on A, A, and E~ production in centrality selected Pb+Pb reactions at 40A and 158A GeV 
and in near-central C+C and Si+Si reactions at 158^1 GeV are presented. Transverse mass spectra, 
rapidity spectra, and multiplicities in dependence of the system size are discussed. Comparisons 
to transport models (UrQMD2.3 and HSD) and to the core-corona approach are performed. While 
A and A production can be described by transport models reasonably well, these models do not 
agree with the H~ measurements. On the other hand, the core-corona picture fits very well the 
system-size dependence of A and H~, while it agrees less well with the A data. 



I. INTRODUCTION 

The production of strange particles has always been 
a key observable in heavy-ion reactions and its enhance- 
ment was one of the first suggested signatures for quark- 
gluon plasma (QGP) formation The predicted en- 
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hancement of strangeness production in nucleus nucleus 
collisions relative to proton-proton reactions was es- 
tablished experimentally some time ago @, Q. It was 
also found that this enhancement is increasing with the 
strangeness content of the particle type 0, [f| . 

However, there are several aspects that make a 
straightforward interpretation of the experimental results 
difficult. One of them is the fact, that the strangeness 
enhancement seems to increase towards lower energies 
[SHU. Another open issue and the topic of this pub- 
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lication is the dependence of the enhancement on the 
system size. The latter is quantified by the number of 
"wounded" nucleons from the colliding nuclei that par- 
ticipate in the reaction. A previous analysis investigated 
the pion, kaon, 0, and A production in (near-) central 
nucleus nucleus collisions of nuclei with different mass 
numbers [§]. It turned out that the enhancement sets 
in already for very small systems and seems to saturate 
for nuclei with mass number above 30, i. e. number of 
participants above 60. The present study extends the 
investigation of the enhancement to non-central Pb+Pb 
collisions and to the production of A and S~ hyperons. 

In the context of statistical models, which generally 
have been quite successful in describing particles yields, 
the experimentally observed strangeness enhancement for 
large systems is due to the increase of the reaction vol- 
ume, which weakens the influence of strangeness conser- 
vation on the production rate [T^. In [ll[ this has been 
modeled by the transition from a canonical ensemble to a 
grand-canonical one. For comparison of the correspond- 
ing statistical model results to experimental data, the 
relation between the size of the ensemble volume V and 
the experimentally accessible number of wounded nucle- 
ons (N w ) has been assumed to be V = (Vb/2) (N w ). 
The parameter Vq, which accounts for the locality of 
the strangeness conservation, is usually fixed to w 7 fm 3 
[Tlj . However, this approach does not provide a satis- 
factory description of the data @. A better agreement 
with global strangeness production at the SPS can be 
achieved by deriving the ensemble volume from a per- 
colation of elementary clusters [T3|. A similar, although 
simplified, line of argument is put forward in the so-called 
core-corona picture [l3l |. Here a heavy-ion collision is 
considered a superposition of a hot and dense core and 
a low density peripheral corona region. While the core 
region corresponds to a large fireball, which experiences 
collective expansion and for which particle production 
should be describable via a large volume canonical en- 
semble, or equivalently by a grand-canonical statistical 
ensemble, the corona is supposed to exhibit the features 
of simple nucleon-nucleon collisions. The contributions 
of the core and the corona depend on both the size of 
the colliding nuclei and the centrality of the collision and 
can be determined via Glauber model calculations [lil ]. 
This concept has recently been applied quite successfully 
to the system-size dependence of strangeness production 
at RHIC energies [11 0J, EE El . 

The comparison of measurements with hadron-string 
transport models, such as UrQMD or HSD, provides in- 
sight into the question whether nucleus-nucleus collisions 
can be described in a hadronic scenario or whether a con- 
tribution from an additional partonic phase is needed. 
Even though these models are not able to describe the 
enhancement of multi-strange particle yields in central 
nucleus-nucleus collisions [8( , their predictions are gener- 
ally close to the measured values for particles containing 
a single strange quark. Their comparison to the mea- 
sured system-size dependence might therefore reveal ef- 
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FIG. 1: (color online) The invariant mass distributions of all 
measured A, A, and H~ candidates in two centrality classes of 
Pb+Pb collisions at 158A GeV. The upper row (a - c) shows 
the most peripheral centrality class, the lower row (d - f ) the 
most central one. The full curves represent a fit to signal 
and background as described in the text, while the dashed 
curves show the background only. The vertical lines denote 
the literature values of the masses 12011. 



fects that go beyond the dominating influence of the re- 
action geometry, such as multi-pion fusion processes and, 
in the case of anti-baryons, absorption in dense nuclear 
matter. 



II. DATA ANALYSIS 

A. Experimental setup and data sets 

The data were taken with the NA49 large acceptance 
hadron spectrometer at the CERN SPS. A detailed de- 
scription of the apparatus can be found in p^ |. With 
this detector tracking is performed by four large-volume 
Time Projection Chambers (TPCs) in a wide range of 
phase space. Two of these are positioned inside two 
superconducting dipole magnets. In order to ensure a 
similar detector acceptance in the center-of-mass system 
for all datasets, the magnetic field was reduced for the 
40A GeV Pb beam by a factor of 4. A measurement of 
the specific energy loss dE/dx in the TPC gas with a 
typical resolution of 4 % provides particle identification 
at forward rapidities. Time-of- flight detectors improve 
the particle identification at mid-rapidity. The centrality 
of a given reaction is determined via the energy measured 
in the projectile fragmentation region by the Zero Degree 
Calorimeter (ZDC) positioned downstream of the TPCs. 
A collimator in front of the ZDC reduces the acceptance 
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FIG. 2: (color online) The pt-integrated reconstruction ef- 
ficiency (a) and feed-down contribution (b) at mid-rapidity 
(A, A: \y\ < 0.4, H~: \y\ < 0.5) as a function of centrality 
as characterized by the measured charged track multiplicity 
for minimum bias Pb+Pb reactions at 40^4 and 158A GeV. 
Open symbols correspond to the values for central Pb+Pb 
collisions The dashed lines represent the A efficiency or 
feed-down contribution respectively, at 158A GeV, resulting 
from a different analysis strategy (see text). 



of the calorimeter to the phase space of the projectile 
fragments and spectator nucleons. 

We present in this paper an analysis of centrality se- 
lected Pb+Pb events taken with a minimum bias trigger 
at beam energies of 40A and 158A GeV in the years 1999 
and 2000, and of near-central C+C and Si+Si events 
measured at 158A GeV in the year 1998. The proper- 
ties of the different datasets are summarized in Tables Q] 
andHTl While for the Pb+Pb collisions the primary SPS 
beam was used, the C and Si ions were created by frag- 
menting the original Pb beam. By tuning the magnetic 
rigidity in the beam line (Z/A = 0.5) and analyzing the 
specific energy loss in scintillation detectors, the corre- 
sponding ion species could be selected. In the following, 
the carbon beam is defined as Z = 6 ions and the silicon 
beam as a mixture of Z = 13 — 15. Two carbon targets 
with thicknesses of 3 mm and 10 mm (561 mg/cm 2 and 
1840 mg/cm 2 , respectively) and a silicon target with a 
thickness of 5 mm (1170 mg/cm 2 ) were used. Further 
details on the analysis of the C+C and Si+Si datasets 
can be found in [9J] . For the study of the minimum bias 
Pb+Pb interactions targets with a thickness of 200 /im 
(224 mg/cm 2 ) were installed. The minimum bias trigger 
is defined by a gas Cherenkov counter that vetoes non- 
interacting projectiles. Centrality classes CO - C4 are de- 
fined by consecutive intervals of spectator energy as mea- 
sured in the ZDC calorimeter. Simulated events from the 
VENUS 4.12 event generator [2l[ were used to relate this 
energy to the number of wounded nucleons (N w ) as given 
by the implemented Glauber model calculation (3] ■ The 
background from non-target interactions is substantially 
reduced by applying cuts on the reconstructed position 
of the primary vertex. After these cuts the centrality 
classes CO - C2 are free of background events, while the 
more peripheral classes C3 and C4 have a contamination 



of less than 2 % and 5 %, respectively. 



TABLE I: Summary of the analyzed Pb+Pb datasets. The 
centrality is quantified by the fraction of the total inelastic 
cross section. (Aw) is the average number of wounded nucle- 
ons per event and a((N w )) the widths of the corresponding 
distributions. For the 158^4 GeV dataset also the fraction 
/((-Aw)) of nucleons that scatter more than once is given [l8| ]. 
AEvcnt is the number of accepted events. 
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TABLE II: Summary of the analyzed near-central C+C and 
Si+Si datasets. The centrality is quantified by the fraction of 
the total inelastic cross section. (N w ) is the average number 
of wounded nucleons per event and Advent the number of 
accepted events. 
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B. A, A, and S reconstruction 

The reconstruction of A, A, and S~ follows the same 
procedures as employed in a previous analysis 0]. A 
detailed description of the methods, together with a list 
of all applied cuts, can thus be found there. Here we just 
summarize the basic principles. 

A and A hyperons were reconstructed from their 
charged decays A — > n~ + p and A — > n + + p (branch- 
ing ratio 63.9 % Candidate pairs were formed by 
combining all reconstructed tracks of positively with all 
tracks of negatively charged particles. Pairs with a dis- 
tance of closest approach (DCA) of less than 0.5 cm any- 
where between the position of the first measured points 
on the tracks and the target plane are treated as possible 
V candidates. The (anti-)protons are identified via their 
specific energy loss (dE/dx) in the TPCs which removes 
the background resulting from a wrong mass assignment. 
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distributions for two centrality classes of Pb+Pb colli- 
sions at 158,4 GeV together with the corresponding fitted 
curves. 



C. Correction for acceptance, reconstruction 
inefficiency and feed-down 

Detailed simulations were made to correct the yields 
for geometrical acceptance and inefficiencies of the re- 
construction procedure. As input to these simulations 
samples of A, A and S - were generated with m t distri- 
butions according to: 



1 dN / m t \ 

m t dmt ^ V T ) 



(1) 



FIG. 3: (color online) The differences between the fully cor- 
rected dN/dy values at mid-rapidity for A at 158A GeV as 
extracted with the standard procedure and the alternative 
analysis strategy (see section III C|l for the different centrality 
bins. The gray area illustrates the systematic error estimate. 



The combinatorial background is further reduced by ap- 
plying additional cuts to V° candidates. These include a 
minimal distance of the reconstructed decay vertex posi- 
tion to the main interaction vertex and the requirement 
that the reconstructed momentum vectors of the A (A) 
candidates should point back to the interaction vertex. 
Due to the lower multiplicities in C+C and Si+Si reac- 
tions, the combinatorial background for A and A is much 
smaller. Therefore, the cuts have been relaxed compared 
to the ones in the analysis of the Pb+Pb data. S~ candi- 
dates were identified via the decay channel EE - — > A + 7r~ 
which has a branching ratio of 99.9 % [20j] . To reconstruct 
the S - , A candidates were selected in an invariant mass 
window of 1.101 - 1.131 GeV/c 2 and combined with all 
measured negatively charged particles in the event. The 
A candidates were subjected to the same selection crite- 
ria as used in the A analysis, except for the momentum 
pointing cut. The reconstructed S~ candidates were re- 
quired to point back to the interaction vertex. In order 
to further reduce the number of fake S - , the trajectories 
of the negatively charged pions from A and S~ decays 
were intersected with the target plane and the distances 
of the resulting positions to the main interaction vertex 
were required to be larger than a minimum value. 

The invariant mass spectra were measured in bins of p t , 
(?7i t — m-o), V, a s well as centrality, and fitted to the sum of 
a polynomial and a signal distribution, the latter deter- 
mined from simulation. The raw yields of A, A, and S~ 
were obtained by subtracting the fitted background and 
integrating the remaining signal distributions in a mass 
window of ±11 MeV/c 2 (±10 MeV/c 2 ) around the nom- 
inal A (S) mass. Typical mass resolutions for A (3), as 
obtained from a fit with a Gaussian, are cr m « 2 MeV/c 2 
(4 MeV/c 2 ). Figure [T] shows examples of invariant mass 



Here pt is the transverse momentum, m t = \/p\ + m 2 , 
and T the inverse slope parameter. In case of the Pb±Pb 
collisions the T parameter was determined by measure- 
ments for central Pb±Pb reactions [H, while for the C±C 
and Si+Si collisions T was set to 200 MeV. The rapidity 
spectra of A and S~ for Pb+Pb reactions were modeled 
by single Gaussian distributions. The flatter A rapidity 
spectra were instead approximated by the sum of two 
(Pb+Pb at 40,4 GeV) or three (Pb+Pb at 158A GeV) 
Gaussians, respectively. For the C+C and Si+Si colli- 
sions the distributions for A and A were assumed to be 
constant as a function of rapidity. The Geant 3.21 pack- 
age [22| was used to track the generated particles and 
their decay products through the NA49 detector. The 
TPC response, which takes into account all known de- 
tector effects, was simulated by software developed for 
the NA49 experiment. The simulated signals were added 
to those of real events on the raw data level and subjected 
to the same reconstruction procedure as the experimen- 
tal data. By determining the fraction of the generated A, 
A and S - which traverse the detector, survive the recon- 
struction, and pass all analysis cuts, the combined accep- 
tances and efficiencies were derived. The corresponding 
correction factors were calculated in bins of pt, (m t — mo), 
y, as well as centrality in the case of Pb+Pb collisions. 

Figure^ summarizes the centrality dependence of the 
efficiencies, including acceptance. While for A and A at 
40A GeV only a 30 % variation can be observed, the dif- 
ference between very peripheral and very central bins is 
a factor of 3 at 158 A GeV. For the A also an alternative 
analysis strategy was employed: on top of the standard 
cuts, only those A were accepted whose daughter tracks 
lie entirely outside the high track density region (CUT-B 
in Q). This reduces the overall efficiency, but has the 
benefit that the track multiplicity effects are slightly re- 
duced compared to the standard analysis (see dashed line 
in Fig. [5fi). The influence of the detector occupancy on 
the efficiency is much smaller for C+C and Si+Si reac- 
tions at 158A GeV. It was found that the reduction of 
the efficiency due to other tracks is ~ 5 % for Si+Si and 
negligible for C+C. Therefore, it was only corrected for 
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FIG. 4: (color online) The transverse mass spectra of A (a), A (b), and H~ (c) at mid-rapidity (A/A: \y\ < 0.4, H~: \y\ < 0.5) 
for Pb+Pb reactions at 40A GeV in different centrality bins. Some of the data points are scaled for clarity. Only statistical 
errors are shown. The solid/dashed lines represent fits with an exponential, where the solid parts denote the mt ranges in 
which the fits were performed. 



in the case of Si+Si reactions. Due to the relaxed analy- 
sis cuts, the efficiencies are generally higher for the small 
systems than for minimum bias Pb+Pb reactions (A/ A: 
- 55 %, 3~: ~ 6 %, both mid-rapidity) [HIH- 

In addition to the directly produced particles, the mea- 
sured yield of A and A contains also contributions from 
the decay of heavier hyperons. The A (A) resulting from 
electromagnetic decays of S° (E°) cannot be separated 
from the directly produced ones via a secondary vertex 
measurement. Thus the presented yields always repre- 
sent the sum A + £° (A + E°). The contribution to A (A) 
from weak decays, however, depends on the chosen anal- 
ysis cuts, because these decay products originate from 
decay vertices with a sizable distance from the main in- 
teraction point. Since the NA49 acceptance for A (A) 
favors those that decay at larger distances, the contribu- 
tion of feed-down A (A) can be quite substantial. There- 
fore, we calculated a correction for the feed-down from 
3~ + S° (S + + S°) decays to the measured A (A) sample 
using the same simulation procedure as described above 
for the efficiency correction. In this case a sample of 
E~ and S° (S + and S°) was generated as input. The 
feed-down correction was then calculated in bins of pt, 
(m t — mo), y, as well as reaction centrality, as the frac- 
tion of reconstructed A (A) which originate from S~ + 5° 
(3+ +3°) decays and pass the same analysis cuts. The 
S~ yields used as input to this procedure are based on 
the measurements presented in this publication, which 
were interpolated to all centrality bins. The extrapola- 
tion of the 5 + yields to the different centrality bins was 
based on the yield measured in central Pb+Pb reactions 



[§[ . It was assumed that the centrality dependence of the 
3 + yields is the same as measured for the 3~. In both 
cases the shapes of the input rapidity and p t distributions 
are a parametrization of the spectra measured for central 
collisions. For the 5° (S°), which are not measured, the 
same shape of the phase space distributions as for the S~ 
was assumed. Their yields are calculated by scaling those 
of S~ (S + ) by the S°/S~ (S°/S + ) ratios taken from sta- 
tistical model fits [25|]. As an example, the pt-hitegrated 
feed-down contribution around mid-rapidity is shown in 
Fig. [2Jd. Since the S yields have a stronger centrality de- 
pendence than the A yields, the feed-down contribution 
changes with centrality. In the case of C+C and Si+Si re- 
actions all yields entering the calculation of the feed-down 
are based on statistical model fits (25|. The parameters 
for their phase space distributions are adjusted such as 
to interpolate between p+p and Pb+Pb reactions. The 
corrections amount to 9 % (10 %) for A and 15 % (20 %) 
for A in C+C (Si+Si) gj. 

D. Systematic errors 

The contributions to the systematic error of the 
AN / dy values measured in centrality selected Pb+Pb re- 
actions are listed in Table IIII1 The first two, the uncer- 
tainty of the background subtraction and the efficiency 
correction, are identical to the ones determined for the 
analysis of the central Pb+Pb datasets @- However, 
there is a difference in the contribution from the feed- 
down correction to the systematic error in the central 
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FIG. 5: (color online) The transverse mass spectra of A (a), A (b), and H~ (c) at mid-rapidity (A/A: \y\ < 0.4, H~: \y\ < 0.5) 
for Pb+Pb reactions at 158A GeV in different centrality bins, and in near-central C+C and Si+Si collisions at 158A GeV. 
Some of the data points are scaled for clarity. Only statistical errors are shown. The solid/dashed lines represent fits with an 
exponential, where the solid parts denote the mt ranges in which the fits were performed. The A spectra for C+C and Si+Si 
are taken from 



TABLE III: Summary of the systematic errors on the 
dN/dy values for minimum bias Pb+Pb reactions. 





■^beam 


Bgnd. 


Eff. 


pt Extra- 


Feed. 


Quad. 




A GeV 


subtr. 


corr. 


polation 


corr. 


sum 


A 


40 
158 


3% 


10% 


3% 


3% 


11% 
11.5% 


A 


40 
158 


3% 


10% 


3% 


8% 


13% 
13.5% 




40/158 


3% 


10% 


3% 




11% 



Pb+Pb analysis. The yields of S~ and particularly of 
S + , arc less precisely measured for non-central Pb+Pb 
reactions than for central ones. While the feed-down con- 
tribution from S~ and S° to A can still be constrained 
reasonably well with the measurement presented here, 
the feed-down estimate from S + and S° to A relies on an 
extrapolation of the measurement in central collisions as- 
suming the same centrality dependence. By varying the 
S input yields inside the errors obtained for the central 
data set and assuming different centrality dependences 
(e.g. scaling with (N w )) the contribution to the system- 
atic error was evaluated. As a result a systematic error 
of 3 % was assigned to the A, while for the A it is 8 %. 
Since the minimum bias data at 158A GeV, and for S - 
also at 40A GeV, do not allow to measure the pt-range 
down to pt = GeV/c, an extrapolation has to be used. 
By comparing the result for the standard fit (exponen- 



tial, as shown in Figs. [4] and [5]) to an extrapolation using 
a fit with a hydrodynamically inspired blast-wave model 
a systematic error of 3 % was determined for these 
cases. Figure 131 demonstrates the consistency of the sys- 
tematic error estimated for the A yield with the spread 
of results obtained from the two analysis strategies dis- 
cussed above. 

The extrapolations in rapidity, which are needed to 
extract the total multiplicities, introduce additional sys- 
tematic errors. The data allow to constrain the widths of 
the fits, as shown in Figs. [7] and [51 only to a certain extent 
which translates into an uncertainty of the extrapolation. 
In the case of A at 158A GeV, the shapes of the rapidity 
distributions are not measured. Therefore a set of as- 
sumptions based on other measurements as discussed in 
[1] was used. The shaded areas in Fig. [5] represent the 
uncertainty of the extrapolations that were included in 
the systematic error. 



III. RESULTS 

A. Transverse mass spectra 

The transverse mass spectra of A, A, and S~ , measured 
around mid-rapidity (A/A: \y\ < 0.4, S~: \y\ < 0.5), are 
shown for different centrality classes of minimum bias 
Pb+Pb collisions at 40A GeV in Fig. Hand at 158A GeV 
in Fig. [5J Also included in Fig. [5] are the mt spectra 
for near-central C+C and Si+Si reactions at 158A GeV. 
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FIG. 6: (color online) The {mt) — mg values at mid-rapidity 
(p/p: -0.5 < y < -0.1, A/A: |y| < 0.4, S - : \y\ < 0.5) for 
Pb+Pb collisions at 40A and 158A GeV, as well as for near- 
central C+C and Si+Si reactions at 158A GeV. The sys- 
tematic errors are represented by the gray boxes. The filled 
symbols correspond to the results obtained from the minimum 
bias data sets, while the open ones represent the (near-)central 
reaction systems. The (anti-)proton data are taken from [27j . 
Also shown are the results from a fit for A and A with the 
core-corona approach (solid lines). 



for the hyperon data compared to previously published 
proton and antiproton results [27| . The mass differences 
between the shown particle species are not very large and 
their (mt) — mo values agree within errors for each par- 
ticular collision system. However, there is a significant 
system-size dependence. A strong increase is observed for 
(N w ) below 60, whereas above this region the values of 
(mt) — mo rise slowly with centrality in Pb+Pb collisions. 



B. Rapidity spectra 

The rapidity spectra of A and A for Pb+Pb collisions 
at 40A and 158A GeV are summarized in Figs. [7] and [5] 
For S~ the statistics of the minimum bias datasets was 
not sufficient to extract rapidity spectra. While the A 
distributions at A0A GeV and the A distributions at 40A 
and 158A GeV have a Gaussian shape at all centralities, 
the A distributions at 158 A GeV are rather flat over the 
measured rapidity range, similar to what has been ob- 
served for central Pb+Pb reactions @. The A and A 
spectra for the near-central C+C and Si+Si collisions are 
shown in Fig. [9] Also here the rapidity distributions for 
A are relatively flat with an indication for a minimum at 
mid-rapidity, which appears to be even more pronounced 
than in the case of the Pb+Pb data. The corresponding 
rapidity densities d N / dy around mid-rapidity for all data 
sets are listed in Tables HVl [VI and ED 



The determination of total multiplicities requires an 
extrapolation into the unmeasured y regions. For this 
purpose the A distributions at 40A GeV were fitted with 
the sum of two Gauss functions of equal width a which 
are displaced symmetrically by s with respect to mid- 
rapidity: 



dN 
dy 



oc exp 



(y - -s) 2 

2a 2 



exp 



(y + sf 

2a 2 



(2) 



The mt spectra were fitted by an exponential as de- 
fined in Eq. |1| in the transverse mass range m t — mo > 
0.2 GeV/c 2 (Pb+Pb data) and m t - m > 0.0 GeV/c 2 
(C+C and Si+Si data). The resulting inverse slope pa- 
rameters T are listed in Tables IIV1 fVl and IVI1 

For a model independent study of the energy de- 
pendence of mt spectra, the average transverse mass 
(m t ) — mo was calculated. To account for the unmea- 
sured m t range two different parametrizations were em- 
ployed to extrapolate: an exponential function (shown in 
Figs. Q] and [5]) and a blast-wave function [2(| (not shown). 
Both provide a good description of the measured data. 
An estimate of the systematic error is derived from the 
differences between the two approaches. The resulting 
values for (m t ) — mo are listed in Tables ITVl fVl and IVII 

Figure O shows the dependence of (m t ) — m on (iV w ) 



At 158A GeV, the data do not allow to determine the 
shape of the A rapidity spectra outside the plateau re- 
gion around mid-rapidity. We therefore use the same as- 
sumptions on the spectral shape that have been applied 
to the central Pb+Pb data at 158 A GeV Q for all cen- 
trality bins in order to extract total multiplicities. This, 
of course, assumes that the widths of the rapidity distri- 
butions do not change substantially with centrality. For 
the A spectra a single Gaussian provides a reasonable fit 
at both beam energies. 

The fitted RMS y values are tabulated in Tables ITVllVl 
and IVII Figure [TO] displays the system-size dependence 
of RMSy. While for A at 40 A GeV an indication for a 
decrease of the widths with increasing centrality can be 
seen, no significant system-size dependence is observed 
for A at both beam energies. 
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FIG. 7: (color online) The rapidity spectra of A and A for Pb+Pb collisions at 40A GeV for the 5 different centrality bins 
CO - C4. The open symbols show data points reflected around mid-rapidity. The systematic errors are represented by the 
gray boxes. Solid lines are fits to the data points, used to extrapolate the measurements in order to extract total yields. Also 
included are calculations with the HSD model [U, OS] (dotted lines) and the UrQMD2.3 model [HHSl (dashed lines). 



C. Particle yields 



Figure [TJJ shows the system-size dependences of the 
rapidity densities dN/dy at mid-rapidity for A, A, and 
E _ , divided by the average number of wounded nucleons 
{N w }. For (N w ) > 60, dN/dy/(N w ) of A is almost in- 
dependent of the system size, while for smaller systems 
at 158A GeV, corresponding to a (N w ) range not cov- 
ered by Pb+Pb collisions, a significant rise with (N w ) is 
observed. However, one should keep in mind that a di- 
rect comparison of near-central C+C and Si+Si reactions 
to peripheral Pb+Pb collisions is complicated by the fact 
that the surface to volume ratio of these reaction systems 
is different. In the case of A, this ratio seems to be inde- 
pendent from the system size even for very small systems. 
The S - , on the other hand, exhibits a weak (N w ) depen- 
dence. For comparison, data of the NA57 collaboration 
are also included in Fig. [TTJ Generally, the dN/dy values 
of NA57 are higher for non-central Pb+Pb collisions, sim- 
ilar to what has already been found for central Pb+Pb 
0. It appears, however, that the discrepancy becomes 
smaller for peripheral collisions. 

The total multiplicities (N) of A and A at 40A and 
158A GeV, as determined from extrapolations of the ra- 
pidity spectra shown in Figs. [7J [5J and [5] normalized 
by (N w ) are summarized in Fig. 1121 A similar picture 
emerges as for the mid-rapidity dN/dy values. In the 
range (iV w ) > 60 the ratio (V)/(A' W ) is independent of 
the system size. Towards smaller system sizes, covered 



only by C+C and Si+Si collisions at 158yl GeV, a signif- 
icant decrease for A is observed, while in case of the A 
this ratio remains constant throughout. 



IV. DISCUSSION 
A. Comparison to transport models 

Transport models allow to study several effects that 
may influence the system size dependence of strange par- 
ticle production, e.g. multi-meson fusion processes, ab- 
sorption of anti-baryons in the fireball, and the evolution 
of the longitudinal distribution of baryon number. 

Multi-meson fusion processes are a possible mechanism 
to enhance the production of anti-baryons [28| and may 
therefore be important for reaching statistical equilib- 
rium yields of multi-strange anti-baryons [23 ] . The HSD 
model [3(| HH Hl| offers the possibility to include these 
fusion processes. The Figures [TTJ and Q2] show compar- 



isons of HSD and UrQMD2.3 [33J, |34|, |35| to the mea- 
sured yields at mid-rapidity and to the total yields. For 
A the predictions of both models are close to the data. 
However, the system-size dependence, especially for the 
total yields, seems to be better described by HSD. The 
multi-meson fusion processes are naturally most impor- 
tant for A. This explains why the spectra and yields for A 
in Pb+Pb predicted by HSD are higher than those from 
UrQMD2.3, which does not feature these processes. HSD 
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FIG. 8: (color online) The rapidity spectra of A and A for Pb+Pb collisions at 158^4 GeV in 5 different centrality bins CO - C4. 
The open symbols show data points reflected around mid-rapidity. The systematic errors are represented by the gray boxes. 
Solid lines are fits to the data points, used to extrapolate the measurements in order to extract total yields. The shaded areas 
in the A spectra sketches the uncertainty due to the different extrapolations. Also included are calculations with the HSD 
model [13, [13, il| (dotted lines) and the UrQMD2.3 model [lllllll (dashed lines). 



gives thus a better description of the measurements. The 
yield of S~ is underestimated by UrQMD2.3 by factors 
of 2 - 3 for all systems. No HSD calculations for the S~ 
are available yet. 

The system-size dependence of anti-baryon yields 
should also be affected by their possible absorption in 
the surrounding dense matter of the fireball. In this case 
one would expect the measured A yield per wounded nu- 
cleon to go down when comparing the small C+C and 
Si+Si systems with central Pb+Pb collisions. In fact, 
the data on {N)/{N W ) for A at 158A GeV seem to ex- 
hibit the expected tendency to decrease from C+C to- 
wards Pb+Pb collisions (see Fig. [T2Ti). quite in contrast 
to the A, where (N)/(N W ) is rather increasing in the re- 
gion (N w ) < 60 (see Fig. [T2b). A similar behavior is 
predicted by UrQMD2.3, where it is, however, stronger 
at AO A GeV than at 158A GeV. But due to the size of the 
systematic error of the measurements, no final conclusion 
can be made whether A production is really affected by 
absorbtion. 

Figures [HJ and [5] include a comparison of the trans- 
port model predictions to the measured rapidity distribu- 
tions of A and A. In the case of A the predicted widths of 
the rapidity distributions from both models, UrQMD2.3 
and HSD, fit the data for all studied systems and energies 
reasonably well (see Fig. [T0|) . A rapidity spectra, which 
are sensitive to the final distribution of baryon number, 
exhibit a significant dependence of their shape on sys- 
tem size. Also here the agreement to the models is fairly 



good at both energies, even though UrQMD2.3 predicts 
a Gaussian shaped distribution at 158A GeV, while the 
data would rather suggest a plateau inside the measured 
region. HSD, on the other hand, describes this flat shape 
relatively well. Similar observations have been made in 
the case of proton rapidity distributions in minimum bias 
Pb+Pb reactions at 158A GeV, where HSD also gives a 
better agreement with the observed flat proton spectra 
than UrQMD2.3 [361 ] ■ The reason for this difference lies 
in a different assumption on when a nucleon is allowed 
to interact again after its first collision. On top of a for- 
mation time of r = 0.8 fm/c, which is implemented in 
both models, HSD requires that the local energy density 
falls below I GeV/fm 3 , which is considered as the criti- 
cal energy density for a phase transition to a QGP. Thus, 
the data would suggest that this additional criterion is 
needed to properly describe the redistribution of baryon 
number in longitudinal phase space due to stopping. 



B. Core-Corona approach 

In order to compare the core-corona approach with the 
data presented here, we generalize the prescription given 
in [l8j and parametrize the system-size dependence of 
any observable X by: 



X((JV W » = (JV W > [/«2Vw» Xcore 

+ (1 - f((N w ))) X c 



(3) 
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FIG. 9: (color online) The rapidity spectra of A and A for 
near-central C+C and Si+Si collisions at 158^4 GeV (the A 
spectra have already been published in Open symbols 
show data points reflected around mid-rapidity, while the sys- 
tematic errors are represented by the gray boxes. Solid lines 
are fits to the data points, used to extrapolate the measure- 
ments in order to extract total yields. Also included are calcu- 
lations with the UrQMD2.3 model [H, H IH| (dashed lines). 
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FIG. 10: (color online) The RMS widths of the rapidity dis- 
tributions RMSy of A and A at 40A and 158A GeV as a 
function of the number of wounded nucleons {N w }. The sys- 
tematic errors are represented by the gray boxes. Filled sym- 
bols correspond to centrality selected data, obtained with a 
minimum bias trigger, while the open ones represent the re- 
action systems measured with a (near-)central trigger. Also 
included are calculations with the HSD model [30l . p3ll . H2] 
(dotted lines) and the UrQMD2.3 model [H, 11 EE[ (dashed 
lines). 



The quantity X can either be the average transverse mass 
(mt) — mo, the rapidity density dN/dy, or the total mul- 
tiplicity (N). The function f({N w )) is here defined as 
the fraction of all participating nucleons, which interact 
more than once, and can therefore be attributed to the 
core region. Since the corona should behave like inde- 
pendent nucleon-nucleon collisions, the quantity A" corona 
corresponds to results of measurements in p+p collisions. 
Thus, the function f((N w )) provides a natural inter- 
polation between p+p and Pb+Pb reactions. We use 
values for f((N w )) (see Table HJ, that have been calcu- 
lated within a Glauber approach for Pb+Pb collisions at 
158 A GeV and have also been used in the toy model com- 
parison discussed in [l8[ . Since the nucleon-nucleon cross 
section changes only slightly between 40 and 158 GeV 
beam energy, we use the same values of f((N w )) for the 
comparison to the 40A GeV data. It should be noted, 
though, that the direct comparison of the curves shown 
here to semi-central C+C and Si+Si collisions is not en- 
tirely correct, since their surface to volume ratio is dif- 
ferent from that in Pb+Pb collisions. This, in princi- 
ple, would require a calculation of f((N w )) specifically 
for these reaction systems. More insight could also be 
gained by studying the smaller systems in several central- 
ity bins, similar to the study of Cu+Cu in [17|. However, 



our available statistics for C+C and Si+Si unfortunately 
does not allow this. 

Based on the above recipe, the system-size dependence 
of (mt) — too for A and A was constructed (solid lines in 
Fig. The (mt) — to values for A in p+p collisions 
are based on an interpolation of p+p data measured at 
various beam energies [38] . The p+p value for A was 
assumed to be the same as for A, since not enough data 
is available to do the extrapolation. The core contribu- 
tions to (mt) — too were adjusted to the measurements 
for central Pb+Pb collisions. In fact, the model provides 
a reasonable description of the measured system-size de- 
pendence in all cases. 

Similarly, the system-size dependence of dN/dy/{N w ) 
and (N) / (N w ) can be predicted using the core-corona 
approach as given by Eq. The solid lines in Figs. [TT1 
and [T121 are based on the same function f({N w )) as has 
been used for (mt) — too. Here, A" corona is adjusted to the 
yields derived from an interpolation of A and A yields 
measured in p+p collisions at different beam energies 
. For the S~ at 158A GeV a preliminary p+p mea- 
surement by NA49 was used 40], while for 40 A GeV no 
p+p input is available so that no comparison to S~ is 
possible at this energy. AT corc is defined in all cases by 
the measured dN/dy, resp. (N), for central Pb+Pb col- 
lisions. The agreement is good for the yields of A and 
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FIG. 11: (color online) The rapidity densities AN /Ay divided 
by the average number of wounded nucleons (N w ) of A, A, 
and E~ at mid-rapidity (A/A: \y\ < 0.4, E~: \y\ < 0.5) for 
Pb+Pb collisions at 40yl and 158A GeV, as well as for near- 
central C+C and Si+Si reactions at 158A GeV, as a function 
of (Nth). The systematic errors are represented by the gray 
boxes. Filled symbols correspond to the minimum bias trig- 
ger, while the open ones represent the online triggered (near- 
central reaction systems. Also shown are data of the NA57 
collaboration [ a, 1371] (open stars) and calculations with the 
HSD model HalSlli] (dotted lines), the UrQMD2.3 model 
[33I.I34II35I (dashed lines), and the core-corona approach (solid 
lines). 



S~ (see: Figs. fTTI and [T2| . However, for A at 158A GeV 
the yields measured in C+C and Si+Si collisions are at 
the same level as for Pb+Pb reactions. This behaviour 
cannot be fitted by the core-corona approach and would 
therefore indicate that also other mechanisms, such as 
A-absorption, need to be taken into account to arrive at 
a proper description of the system-size dependence. 



V. SUMMARY 

A measurement of A, A, and S~ production in central- 
ity selected Pb+Pb collisions at 40A and 158A GeV and 
in near-central C+C and Si+Si collisions at 158A GeV 
is presented. The first moments of the transverse mass 
spectra ((mt) — too) exhibit only a weak system-size de- 
pendence for (N w ) > 60, while for the small systems a 
rapid rise of (m t ) — too with increasing system size is 
observed. The rapidity distributions of A at 40^4 GeV 
and of A at 40A and 158^4 GeV have a Gaussian shape. 
For A at 158 A GeV the rapidity spectra are rather flat 
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FIG. 12: (color online) The total multiplicities (N) divided by 
the average number of wounded nucleons {N w } of A and A for 
Pb+Pb collisions at 40A and 158^1 GeV, as well as for near- 
central C+C and Si+Si reactions at 158A GeV, as a func- 
tion of {N w ). The systematic errors are represented by the 
gray boxes. Filled symbols correspond to the minimum bias 
trigger, while the open ones represent online triggered (near- 
Central reaction systems. Also shown are calculations with 
the HSD model [p, [U, El (dotted lines), the UrQMD2.3 
model [33|, I34I |35|] (dashed lines), and the core-corona ap- 
proach (solid lines). 



in the measured region —1.6 < y < 1.2. Generally, 
no pronounced system-size dependence of the widths 
of the rapidity distributions is observed. Only the A 
spectra at 40A GeV might show some indication for a 
slight narrowing with increasing centrality. The mea- 
sured AN I dy j (N w ) and (N)/{N W ) values rise rapidly 
with system size for small systems ((N w ) < 60) and do 
not change very much any more from then on. The core- 
corona approach describes the system-size dependence 
of (mt) — too for all particle species discussed here. It 
also reproduces the system-size dependence of the mid- 
rapidity dN/dy values and of the total multiplicities of A 
and S~ . However, the A spectra measured for C+C and 
Si+Si collisions at 158A GeV suggest a flatter system- 
size dependence than expected in the core-corona picture. 
Generally, the results of the hadronic transport models 
UrQMD2.3 and HSD for A and A are close to the data, 
with the exception of the underprediction of the abso- 
lute A yields at both energies by UrQMD2.3. But both 
models predict a system-size dependence of the A total 
multiplicity similar to the measurement. This might indi- 
cate that absorption of A in the dense hadronic medium, 
which is taken into account in the hadronic transport 
models, has a visible effect. However, these models are 
not able to describe the production of baryons with mul- 
tiple strangeness. UrQMD2.3, for example, underesti- 
mates the yields of S~ to a large extent (factor 2-3). 
The HSD model seems to provide a better description 
of the A rapidity spectra than UrQMD2.3 due to an im- 
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proved implementation of the stopping mechanism. 
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TABLE IV: The rapidity densities at mid-rapidity (A/A: \y\ < 0.4, H": \y\ < 0.5), the total multiplicities (N), the RMS 
widths of the rapidity distributions RMSy calculated from the fits shown in Fig. [7] the average transverse masses (rat) — mo, 
and the inverse slope parameters T for Pb+Pb collisions at 40A GeV. The first error is statistical, the second systematic. 





Centrality 
class 


(iVw) 


dN/dy 


(N) 


RMSy 


(mt) - m 
(MeV/c 2 ) 


T 

(MeV) 


A 





351+3 


14.9+0.3+1.6 


41.1±0.8±5.3 


0.96±0.06±0.11 


324±10±33 


268± 7±16 




1 


290+4 


12.6+0.2+1.4 


34.1±0.6±4.4 


0.98±0.07±0.11 


310± 8±31 


256± 4±15 




2 


210+6 


8.4+0.1+0.9 


23.6±0.3±3.1 


0.99±0.05±0.11 


299± 6±30 


247± 3±15 




3 


142+8 


5.3+0.1+0.6 


16.4±0.2±2.1 


1.14±0.13±0.13 


277± 6±28 


230± 4±14 




4 


93+7 


3.2+0.05+0.4 


9.9±0.1±1.3 


1.12±0.08±0.12 


261± 6±26 


220± 5±13 


A 





351+3 


0.29+0.04+0.04 


0.59±0.08±0.12 


0.79±0.16±0.10 


405±72±53 


325±81±32 




1 


290+4 


0.23+0.03+0.03 


0.43±0.05±0.09 


0.71±0.12±0.09 


366±56±48 


299±40±30 




2 


210+6 


0.18+0.02+0.02 


0.35±0.03±0.07 


0.69±0.09±0.09 


321±42±42 


276±38±28 




3 


142+8 


0.14+0.01+0.02 


0.27±0.02±0.05 


0.73±0.08±0.10 


392±36±51 


346±69±35 




4 


93+7 


0.09+0.01+0.01 


0.17±0.01±0.03 


0.77±0.10±0.10 








0-1 


314+4 


1.25+0.23+0.14 






286±33±34 


232±19±14 




2-3 


178+8 


0.49+0.07+0.05 






272±23±33 


233±18±14 


TABLE V: The rapidity densities at mid-rapidity (A/A: \y\ < 0.4, H~: \y\ < 0.5), the total multiplicities (N), the RMS widths 
of the rapidity distributions RMSy calculated from the fits shown in Fig. [8] the average transverse masses (mt) — tuq, and the 
inverse slope parameters T for Pb+Pb collisions at 158A GeV. The first error is statistical, the second systematic. 




Centrality 
class 


(iVw) 


dN/dy 


(N) 


RMSy 


(mt) - m 
(MeV/c 2 ) 


T 

(MeV) 


A 





352+3 


12.9+0.7+1.5 


55.3±1.8±8.8 




402±43±48 


346±34±21 




1 


281+4 


10.0+0.4+1.2 


45.9±1.0±7.3 




354±21±43 


296±14±18 




2 


196+6 


6.9+0.2+0.8 


32.0±0.5±5.1 




361±16±43 


307±11±18 




3 


128+8 


4.4+0.1+0.5 


19.7±0.3±3.2 




353±15±42 


303±12±18 




4 


85+7 


2.7+0.1+0.3 


11.9±0.2±1.9 




316±14±38 


274±15±16 


A 





352+3 


1.4+0.3+0.2 


3.6±0.4±0.6 


0.96±0.14±0.13 


580±148±75 


507±211±51 




1 


281+4 


0.92+0.14+0.12 


3.0±0.3±0.5 


1.11±0.15±0.15 


443±109±58 


372±62±37 




2 


196+6 


0.78+0.07+0.11 


2.2±0.1±0.4 


0.98±0.08±0.13 


345±31±45 


296±27±30 




3 


128+8 


0.54+0.04+0.07 


1.4±0.1±0.2 


0.99±0.10±0.13 


345±27±45 


302±34±30 




4 


85+7 


0.36+0.03+0.05 


0.8±0.05±0.14 


0.83±0.07±0.11 


340±23±44 


309±50±31 




0-1 


309+4 


1.43+0.33+0.16 






317±39±38 


244±41±15 




2-3 


164+8 


0.59+0.10+0.06 






327±29±39 


264±39±16 




4-5 


54+7 


0.17+0.03+0.02 






333±29±40 


261±37±16 




0-2 


262+4 


0.96+0.10+0.11 






330±24±40 


263±14±16 


TABLE VI: The rapidity densities at mid-rapidity (A/A: \y\ < 0.4, H": \y\ < 0.5), the total multiplicities (N), the RMS 
widths of the rapidity distributions RMSy calculated from the fits shown in Fig. [9] the average transverse masses (mt) — mo, 
and the inverse slope parameters T for near-central C+C and Si+Si collisions at 158A GeV. The first error is statistical, the 
second systematic. 




Reaction 
system 


(N w ) 


dN/dy 


(N) 


RMSy 


(mt) - m 
(MeV/c 2 ) 


T 

(MeV) 


A 


C+C 


14+2 


0.24 ±0.01+0.04 


1.32±0.05±0.32 




224± 6±27 


199± 8±15 




Si+Si 


37+3 


0.88 ±0.04+0.13 


3.88±0.16±0.56 




267±16±32 


235± 9±16 


A 


C+C 


14+2 


0.064+0.003+0.010 


0.18±0.02±0.03 


0.92±0.04±0.16 


204± 9±27 


184±11±17 




Si+Si 


37+3 


0.16 ±0.007±0.038 


0.49±0.05±0.11 


0.96±0.13±0.21 


230±10±30 


205± 9±17 




Si+Si 


37+3 


0.07 ±0.01 ±0.01 






239±16±29 


210±13±13 



